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Abstract 

Heat conduction in graphite has been studied for decades because of its exceptionally large 
thermal anisotropy. While the bulk thermal conductivities along the in-plane and cross-plane 
directions are well known, less understood are the microscopic properties of the thermal phonons 
responsible for heat conduction. In particular, recent experimental and computational works 
indicate that the average phonon mean free path (MFP) along the c-axis is considerably larger 
than that estimated by kinetic theory, but the distribution of MFPs remains unknown. Here, 
we report the first quantitative measurements of c-axis phonon MFP spectra in graphite at a 
variety of temperatures using time-domain thermoreflectance measurements of graphite flakes with 
variable thickness. Our results indicate that c-axis phonon MFPs have values of a few hundred 
nanometers at room temperature and a much narrower distribution than in isotropic crystals. 
At low temperatures, phonon scattering is dominated by grain boundaries separating crystalline 
regions of different rotational orientation. Our study provides important new insights into heat 
transport and phonon scattering mechanisms in graphite and other anisotropic van der Waals solids. 
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Thermal transport in anisotropic layered materials has attracted broad interest in 
fundamental science and applications [THT^ due to their large thermal anisotropy. The best 
known of these materials, graphite, possesses ultrahigh in-plane thermal conductivity of 
~ 2000 W/mK[T3] and cross-plane (c-axis) thermal conductivity around 6.8 W/mK [T^ 
at room temperature, yielding an anisotropy factor of around 300. During the past 
decade, research has mainly focused on thermal transport properties along the in-plane 
direction, especially in graphene, the ideal two-dimensional system isolated from the 
ab-plane [T3H22] • Such a high in-plane thermal conductivity has important applications in 
heat dissipation^ ES]- In addition, the low thermal contact resistance of ultra-thin graphite 
and graphene flakes make them promising candidates for thermal interface materials |24]. 

Compared with works on the in-plane direction, thermal transport along the c-axis of 
graphite has not yet been extensively studied. Early works used kinetic theory to estimate 
that the average MFP along the c-axis is only a few nanometers at room temperature ^31EB] . 
Recently, several studies have indicated that the average phonon mean free path along 
c-axis (c-MFP) is actually much higher than this simple estimation. Sadeghia et ah 
estimated the average c-MFP is about 20 nm[27]. Harb et ah reported a cross-plane 
thermal conductivity of only 0.7 W/mK in a 35 nm thick graphite film, giving strong 
evidence of boundary scattering of phonons with MFPs much longer than 35 nm|28]. Using 
molecular dynamics simulations, Wei et alj^H]- observed the presence of very long c-MFP 
phonons on the order of several hundred nanometers. However, this study was limited in 
the thickest samples that could be studied by computational limitations. Experimental 
measurements by Yang et al. suggested that the average MFP along c-axis should be much 
greater than 92 nm at room temperature [11], but this estimation is based on a model of an 
isotropic van der Waals (vdW) solid rather than the actual anisotropic phonon dispersion |llj. 

A pioneering experimental work by Fu et al. has recently provided strong evidence that 
the average c-MFP at room temperature is around 204 nm[T2|. In this work, thickness 
dependent thermal conductivity of graphite ribbons was obtained with the 3a; method and 
the c-MFP estimated by fitting the experimental data[T2|. However, due to experimental 
limitations important questions still remain regarding the spectrum of c-MFPs over a range 
of temperatures. For example, the limited temperature range of the measurements of Fu et 
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al makes identifying the dominant scattering mechanisms in different temperature regimes 
challenging. 

In this work, we report the hrst quantitative measurements of the temperature-dependent 
c-MFP spectrum of graphite. Using a theoretical approach to interpret systematic mea¬ 
surements of c-axis thermal conductivity over a wide range of sample thicknesses and 
temperatures, we reconstruct the c-MFP spectrum at each temperature from which we can 
identify the key scattering mechanisms. We hnd that at room temperature, the c-MFPs 
range from 40 nm to 250 nm, a much narrower distribution than in isotropic crystalline 
materials but far larger than early estimates of a few nanometers. At low temperatures, 
the thermal conductivity is limited by phonon scattering at grain boundaries separating 
crystalline regions of different rotational orientations, yielding MFPs between 100-600 nm. 
Our study provides important new insights into the microscopic processes governing heat 
conduction in graphite. 

We conducted thermal conductivity measurements of thin dims of graphite using a 
thermal characterization technique, time-domain thermoreflectance (TDTR). TDTR is a 
standard characterization method in the thermal sciences that operates by impulsively 
heating a sample with a laser pulse and observing the transient thermal decay with a 
probe beam. The thermal conductivity is obtained by fitting the thermal decay curve to a 
thermal model mEU]. In this study, we use a two-tint implementation of this method as 
described in Ref. EH Briefly, a Ti: sapphire laser oscillator generates a series of femtosecond 
pulses at a repetition frequency of 76 MHz and wavelength of 785 nm. The pump pulse 
impulsively heats the sample, and the time-delayed probe pulse measures the change 
in optical reflectance of the Al transducer due to the temperature change, yielding a 
temperature decay curve. The pump pulse train is modulated at a frequency between 0.5 
MHz to 14 MHz to enable lock-in detection. The laser spot sizes of pump and probe ranged 
from 13.5 pm to 60 pm and 10 pm (hxed), respectively. The measurements were found 
to be independent on pump size. Low temperature measurements were performed in an 
optical cryostat under vacuum of around 10“® torr. 

The samples in this study consisted of a sandwich structure as shown in Figure la, with 
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the A1 transducer film and sapphire substrate sandwiching graphite samples of various 
thicknesses. Graphite flakes were first mechanically exfoliated from bulk HOPG samples 
and placed on top of sapphire wafers, with adhesion being by the vdW force. Large, flat 
regions in graphite samples were carefully chosen via optical microscopy. Their thicknesses 
were characterized by atomic force microscope (AFM) (Supplementary Information). Then, 
the A1 transducer layer was deposited on top of the graphite flakes by electron beam 
evaporation. The thicknesses of the A1 Aims were measured by AFM and had typical 
values of 70 nm to 86 nm. Figure lb shows an optical image of one typical sample in our 
experiment (see additional images in the SI). Samples were carefully selected from a large 
range of thicknesses from 90 nm to 500 pm, with the thickest sample corresponding to 
bulk HOPG. Every sample was measured at different temperatures ranging from 40 K to 
294 K. Gompared with a previously reported process for fabricating similar graphite thin 
film samples, in which a series of wet lithography and dry etching steps were used[T2], our 
approach has several advantages. First, we can fabricate very thick samples (on the order of 
tens of microns) onto the substrates. In contrast, firmly attaching a very thick but narrow 
graphite strip onto the substrate and maintaining this state during lithographic patterning 
can be challenging. In this work, we select large area (on the order of millimeters), but 
generally thin, graphite flakes which can be easily placed on the substrates. Second, by 
avoiding the etching operation, our method eliminates potential damage to the samples 
when their lateral geometries are trimmed. 

Thermal properties of graphite flakes were obtained by fitting measured temperature 
decay curve to a thermal model[7l |30]. Figures Ic-lf show a typical experimental data set 
from our TDTR experiment and the corresponding fitting results from a thermal model, 
demonstrating excellent agreement. The fitting process for such a multi-layer, anisotropic 
structure is more complicated than for traditional TDTR because there are 4, rather than 
2, unknown fitting parameters: the cross-plane thermal conductivity of the graphite flakes, 
denoted as kc] interface thermal conductance between the A1 film and the graphite flake, Gi] 
the anisotropy of the graphite flake, S] and the interface thermal conductance between the 
graphite flake and the sapphire substrate, G 2 . Further, every different thickness of graphite 
is a separate sample that may not necessarily have the same value of G 2 or anisotropy S. 
To overcome these issues, we performed a multi-step fitting process. We first created a 


4 


matrix of possible S and G 2 values based on general knowledge of the range of these two 
parameters. Each matrix element consists of an S and a G 2 value. Then, we £t the data to 
the thermal model by hxing S and G 2 values from each matrix element and varying kf. and 
Gi as the htting parameters. After performing this htting operation for the whole matrix, 
a screening was performed by examining the quality of the fit according to a threshold 
criterion. We recorded the thermal conductivity acquired from each high-quality £t and 
discarded those from low quality fits. The remaining htting results are used for calculating 
the average value and the error bar of the measurement. An example of this procedure is 
shown in the supplementary information. Following this procedure for each sample and 
temperature, we extract the thermal conductivity of the graphite hake versus thickness and 
temperature. 

Figures 2a and 2b show the results for samples thicker and thinner than 400 nm, 
respectively. For samples thicker than 400 nm. Figure 2a, the thermal conductivity is 
largely independent of thickness and exhibits a peak at around 100 K, a common trend 
that is characteristic of intrinsic phonon-phonon scattering at room temperature and 
boundary scattering at low temperature. Considering the measured thermal conductivity is 
independent of thickness, these boundaries must consist of internal grain boundaries. The 
trend of the curve as well as the peak position are consistent with prior experimental and 
theoretical works [221 ES]- The measured bulk thermal conductivity at room temperature 
is between 6.4 W/mK to 7.0 W/mK, in agreement with the accepted value of 6.8 W/mK[TT]. 

Samples with thicknesses less than 400 nm. Figure 2b, have a diherent trend with 
temperature, slowly increasing rather than exhibiting a peak. This trend is characteristic 
of samples that have a thermal conductivity limited by boundary scattering. Because the 
thermal conductivity decreases with decreasing sample thickness, the relevant boundary 
for these samples is the sample boundary, thereby indicating that c-MFPs are on the order 
of the sample thickness. The crossover from the bulk trend to the boundary scattering 
trend occurs at thickness ~ 200 nm. Our observations are consistent with Fu et al’s recent 
experimental work [12]. as shown in Figures 2a and 2b. 

Figures 2c and 2d show the thermal conductivity versus thickness at a series of diherent 
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temperatures. Owing to the wide range of sample thicknesses, we are able to directly 
determine the saturation thickness, dehned as the thickness above which the sample recovers 
its bulk thermal conductivity. One can notice that the saturation thickness increases with 
decreasing temperature, indicating that c-MFPs increase with decreasing temperature, as 
expected. At temperatures less than 100 K, the saturation thickness reaches around 600 nm 
and is independent of temperature. This value is much shorter than the thicknesses of the 
thick samples in this study, which range from microns to hundreds of microns. Therefore, 
dominant scattering mechanism must be intrinsic scattering by defects or grain boundaries 
along c-axis, separately conhrming the original conclusion obtained from the temperature 
dependence of the thermal conductivity. 

Figure 2c also illustrates a comparison with experimental data from Fu et al in a similar 
range of thicknesses at room temperature |12] . showing a quite consistent trend. However, 
due to the limitation of the thickness of the samples used by Fu et al, the maximum of which 
was 714 nm, the saturation value was only estimated via theoretical extrapolation. Here, 
we are able to directly measure the saturation thickness and the bulk value. We notice that 
our measured bulk values are higher than those estimated by Fu et al. One possible reason 
for this difference is that the etching process required for the Sen method may introduce 
additional defects that scatter phonons, lowering the bulk thermal conductivity. Therefore, 
maintaining large sample sizes in the basal plan is still important even for measuring 
cross-plane thermal properties. 

Our measurements show that scattering from internal grain boundaries largely sets the 
c-axis thermal conductivity of graphite at low temperatures. To identify the structures 
responsible for grain boundary scattering, we performed transmission electron microscopy 
(TEM) studies on our samples. Figure 3a, a bright held (BF) TEM image, clearly shows 
the existence of some grains of different contrast from the bulk (dark regions), which 
are separated by a few hundred nanometers on average. Overall, the sample is highly 
crystalline as conhrmed by diffraction patterns of the sample (see SI). These observations 
are consistent with Park et al’s work|34]. 

Figure 3b is a dark held TEM image of the same region in Figure 3a. All the bright 
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regions indicate grains of a specific crystalline orientations, which are distinct from the 
other areas of the sample. These bright regions can be readily located as corresponding 
dark regions in the same areas of Figure 3a, which means that the variations of contrast in 
BF images originates from grains of different crystalline orientations. 

Figure 3c is a magnihed BF TEM image on the same sample. Considerable contrast is 
observed even within the dark regions, indicating differences in crystal structure between 
the bright background and these regions. High resolution TEM images on the boundary of 
dark and bright regions. Figure 3d, shows that both areas demonstrate good crystallinity 
and atomic registry at the interface. These observations imply that the grains maintain the 
same c-axis with the rest of the sample but are of different rotational orientations. Because 
the saturation thickness is on the order of hundreds of nanometers, precisely the spacing of 
these rotational mismatches, our results indicate that the rotational orientation differences 
are the structural feature responsible for phonon scattering. 

With these data, we can qualitatively estimate that the typical phonon c-MFP on the 
order of hundreds of nanometers at room temperature. However, to obtain quantitative 
information on the c-MFP spectrum, we need a way to link the measured thermal 
conductivities to the c-MFP spectrum. We have recently reported a theoretical approach 
that makes this hnk |H3] - IH7] . which we now apply here. 

The cross-plane thermal conductivity reduction in graphite thin Elms can be attributed 
to boundary scattering of phonons [37] in a conceptually similar manner to in-plane thermal 
conductivity reduction described by Fuchs-Sondheimer theory [381139]. In a recent work, we 
derived the equation that relates the cross-plane thermal conductivity to the him thickness 
for an isotropic crystal: 

poo poo 

h= / ^(Kn,)/(A^)dA^ = / L-^K{Kn^)F{A^)dA^ (1) 

Jo Jo 

where Kn^^ = A ^/L is the Knudsen Number, A^ is the c-MFP, L denotes the sample thickness 
along the c-axis, /(A^^) and E(A^) are differential and accumulative c-MFP spectra related 
by E(A^) = f{A)dA, and S'(Knt^) is the heat hux suppression function that describes 
the reduction in phonon MFP due to boundary scattering. The kernel iF(Knt^) is dehned 
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as K{Kn.^) = —(iS'(Kn(^)/(iKn^. The suppression function was originally derived in Ref. |37] 
for an isotropic crystal. Repeating the derivation for a crystal with an arbitrary dispersion 
relation, we hnd 

R(Kn,,k) = 1 - Kne,k(l - e“^) (2) 

where Knc^k = is the Knudsen number along the c-axis, dehned as the component of 

the MFP along the c-axis Ac,k normalized by the thickness of the sample L. Both ^(Knc^k) 
and iF(Knc,k) are plotted in the inset of Figure 4a. 

Using this function, the c-MFP spectra can be recovered from the data by solving 
an inverse problem as described in Ref. 133 The original experimental data and the 
reconstructed c-MFP spectra at different temperatures are shown in Figure 4. At room 
temperature, the c-MFP spectrum has a minimum and maximum of around 40 nm and 250 
nm, respectively. The result is consistent with Yang et al’s qualitative estimation that the 
average of c-MFPs can be well over 100 nm[TT]. As temperature decreases, both bounds 
increase and the bandwidth increases monotonically. At 40 K, the lower and upper bounds 
have increased to about 100 nm and 600 nm, respectively. 

Note that the width of the spectrum is considerably narrower than in common isotropic 
crystals such as SigQl E], which has a MFP spectrum spanning at least four orders of 
magnitude. We attribute the narrow bandwidth to two main factors. First, along the 
c-axis, graphite is a vdW solid that does not support high frequency phonons. As a result, 
the range of phonon frequencies that contribute to c-axis heat conduction is considerably 
narrower than in isotropic crystals. Only phonons whose frequencies are lower than 4 THz 
can effectively transport along the c-axis|32]. The lack of contribution from high frequency 
phonons imposes a restriction on the minimum possible c-MFP. Second, grain boundaries 
reduce the maximum c-MFPs to a value comparable to the grain size. Therefore, the 
c-MFP spectrum is restricted both at small and large c-MFPs. 

We plot the differential c-MFP spectra at various temperatures in Figure 5a, allowing 
the change of the bandwidth and peak of the spectrum with temperature to be clearly 
observed. As temperature decreases, both the upper and lower bounds of the c-MFPs 
increase. Meanwhile, a small broadening of the bandwidth can be observed. 



To more compactly describe these distributions, we define c-MFPp as the c-MFP at 
the maximum of the differential distribution. Additionally, we calculate the midpoint 
MFP, c-MFPmid, which indicates the 50% point of the accumulative distributions in Fig 
4. The two parameters, together with lower and upper c-MFP bounds as functions of 
temperature, are shown in Figure 5b. Although all the four parameters increase with 
decreasing temperature, the upper bound of c-MFP increases faster than the rest. The 
c-MFPp is almost at the first trisection point of between the upper and lower bounds. 
These results provide a clear indication of the average c-MFP and the bandwidth at a range 
of temperatures in graphite. 

Our study of c-MFP spectra of graphite provides important scientific insights as well 
as guidance for engineering graphitic thermal materials for applications such as graphite 
nanoplatelet papers |13l El] and thermal interface materials [IS]. In particular, our work 
shows that the c-axis thermal conductivity is only achieved in graphite films larger 
than 400 nm at room temperature; thinner films will conduct heat less effectively along 
the c-axis than would be expected according to the bulk thermal conductivity due to 
boundary scattering. On the other hand, decreasing the thermal conductivity appears to 
be quite achievable due to the long c-MFPs. One potential approach to achieve lower 
thermal conductivity is introducing nanoparticles in the bulk graphite between layers or 
grains[ini HZ], or, as discussed in this work, introducing additional rotational mismatches 
of grains to scatter phonons. 

Scientifically, our work demonstrates that the microscopic properties of thermal phonons 
can be considerably different in highly anisotropic materials than in isotropic solids. 
In particular, the c-axis MFP spectrum of graphite is much narrower than in isotropic 
crystals due to restrictions on both the maximum and minimum c-MFPs. Additionally, our 
approach of combining thermal measurements over length scales comparable to MFPs with 
TEM characterization shows considerable promise for studying the microscopic processes 
governing heat conduction in other anisotropic solids. 

In conclusion, we have reported the first quantitative measurements of the c-axis MFPs 
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of thermal phonons in graphite. We hnd that at room temperature, the c-MFPs are on 
the order of 100-200 nm and are primarily limited by intrinsic phonon-phonon scattering, 
while at low temperatures boundary scattering from internal grains governs the thermal 
transport. We also hnd that the c-MFP spectrum of graphite is unusually narrow compared 
with that of common isotropic crystals due to the narrow range of phonon frequencies 
that are supported along the c-axis. Our work provides important insights necessary to 
understand and engineer heat conduction in graphite and other vdW solids. 


Acknowledgements 

The authors are grateful to the Kavli Nanoscience Institute at Caltech for the availability 
of critical cleanroom facilities and to the Lewis Group at Caltech for use of certain facilities. 
The authors also thank Matthew H. Sullivan for FIB assistance and Carol M. Garland for 
TEM assistance. This work was supported by a start-up fund from the California Institute of 
Technology and by the National Science Foundation under CAREER Grant CBET 1254213. 

Additional Information 

The authors declare no competing hnancial interests. 


10 



[1] Minnich, A. J. Phys. Rev. B 2015 , 91, 085206. 

[2] Jang, W.; Bao, W.; Jing, L.; Lau, C. N.; Dames, C. Applied Physics Letters 2013 , 103, 133102. 

[3] Dresselhaus, M. MRS Bulletin 1987 , 12, 24-28. 

[4] Chiritescu, C.; Cahill, D. G.; Nguyen, N.; Johnson, D.; Bodapati, A.; Keblinski, P.; Zschack, P. 
Science 2007 , 315, 351-353. 

[5] Yan, Z.; Liu, G.; Khan, J. M.; Balandin, A. A. Nat Commun 2012 , 3, 827. 

[6] Losego, M. D.; Blitz, I. P.; Vaia, R. A.; Gahill, D. G.; Braun, P. V. Nano Letters 2013, 13, 
2215-2219, PMID: 23594105. 

[7] Schmidt, A. J.; Chen, X.; Chen, G. Review of Scientific Instruments 2008, 79, 114902. 

[8] Luckyanova, M. N.; Johnson, J. A.; Maznev, A. A.; Garg, J.; Jandl, A.; Bulsara, M. T.; 
Fitzgerald, E. A.; Nelson, K. A.; Chen, G. Nano Letters 2013 , 13, 3973-3977, PMID: 
23952943. 

[9] Minnich, A. J. arXiv.org 2015, DOI: http://arxiv.org/abs/1508.00611, 

[10] Foley, B. M.; Herndez, S. C.; Duda, J. C.; Robinson, J. T.; Walton, S. G.; Hopkins, P. E. 
Nano Letters 2015 , 15, 4876-4882, PMID: 26125524. 

[11] Yang, J.; Shen, M.; Yang, Y.; Evans, W. J.; Wei, Z.; Chen, W.; Zinn, A. A.; Chen, Y.; 
Prasher, R.; Xu, T. T.; Keblinski, P.; Li, D. Phys. Rev. Lett. 2014 , 112, 205901. 

[12] Fu, Q.; Yang, J.; Chen, Y.; Li, D.; Xu, D. Applied Physics Letters 2015 , 106, 031905. 

[13] Touloukian, Y. Thermophysical Properties of Matter: Thermal conductivity: nonmetallic 
solids, by Y. S. Touloukian and others] The TPRC data series; IFI/Plenum, 1970. 

[14] Taylor, R. Philosophical Magazine 1966 , 13, 157-166. 

[15] Jang, W.; Chen, Z.; Bao, W.; Lau, C. N.; Dames, C. Nano Letters 2010 , 10, 3909-3913, 
PMID: 20836537. 

[16] Ghosh, S.; Bao, W.; Nika, D. L.; Subrina, S.; Pokatilov, E. P.; Lau, C. N.; Balandin, A. A. 
Nat Mater 2010 , 9, 555-558. 

[17] Seol, J. H.; Jo, L; Moore, A. L.; Lindsay, L.; Aitken, Z. H.; Pettes, M. T.; Li, X.; Yao, Z.; 
Huang, R.; Broido, D.; Mingo, N.; Ruoff, R. S.; Shi, L. Science 2010, 328, 213-216. 

[18] Bonini, N.; Garg, J.; Marzari, N. Nano Letters 2012 , 12, 2673-2678, PMID: 22591411. 

[19] Lindsay, L.; Li, W.; Garrete, J.; Mingo, N.; Broido, D. A.; Reinecke, T. L. Phys. Rev. B 2014 , 


11 



89, 155426. 

[20] Paulatto, L.; Mauri, F.; Lazzeri, M. Phys. Rev. B 2013, 81, 214303. 

[21] Bae, M.-H.; Li, Z.; Aksamija, Z.; Martin, P. N.; Xiong, F.; Ong, Z.-Y.; Knezevic, I.; Pop, E. 

Nat Commun 2013, 1734. 

[22] Xu, X.; Pereira, L. F. C.; Wang, Y.; Wu, J.; Zhang, K.; Zhao, X.; Bae, S.; Tinh Bui, C.; 
Xie, R.; Thong, J. T. L.; Hong, B. H.; Loh, K. P.; Donadio, D.; Li, B.; zyilmaz, B. Nat 
Commun 2014, 5, 4689. 

[23] Ghosh, S.; Calizo, L; Teweldebrhan, D.; Pokatilov, E. P.; Nika, D. L.; Balandin, A. A.; Bao, W.; 
Miao, E.; Lau, C. N. Applied Physics Letters 2008, 92, 151911. 

[24] Chen, Z.; Jang, W.; Bao, W.; Lau, C. N.; Dames, C. Applied Physics Letters 2009, 95, 161910. 

[25] Shen, M.; Schelling, P. K.; Keblinski, P. Phys. Rev. B 2013, 88, 045444. 

[26] Tanaka, T.; Suzuki, H. Carbon 1972, 10, 253 ~ 257. 

[27] Sadeghi, M. M.; Jo, L; Shi, L. Proceedings of the National Academy of Sciences 2013, 110, 
16321-16326. 

[28] Harb, M.; von Korff Schmising, C.; Enquist, H.; Jurgilaitis, A.; Maximov, L; Shvets, P. V.; 
Obraztsov, A. N.; Khakhulin, D.; Wulff, M.; Larsson, J. Applied Physics Letters 2012, 101, 
233108. 

[29] Wei, Z.; Yang, J.; Chen, W.; Bi, K.; Li, D.; Chen, Y. Applied Physics Letters 2014, 104, 
081903. 

[30] Cahill, D. G. Review of Scientific Instruments 2004, 75, 5119-5122. 

[31] Kang, K.; Koh, Y. K.; Chiritescu, C.; Zheng, X.; Cahill, D. G. Review of Scientific Instruments 
2008, 79, 114901. 

[32] Ho, C. Y.; Powell, R. W.; Liley, P. E. Journal of Physical and Chemical Reference Data 1972, 
1, 279-421. 

[33] Fugallo, G.; Cepellotti, A.; Paulatto, L.; Lazzeri, M.; Marzari, N.; Mauri, F. Nano Letters 
2014, 14, 6109-6114, PMID: 25343716. 

[34] Park, S.; Floresca, H. C.; Suh, Y.; Kim, M. J. Carbon 2010, 4S, 797 - 804. 

[35] Minnich, A. J. Phys. Rev. Lett. 2012, 109, 205901. 

[36] Zhang, H.; Hua, C.; Ding, D.; Minnich, A. J. Sci. Rep. 2015, 5, 9121. 

[37] Hua, C.; Minnich, A. J. Journal of Applied Physics 2015, 117, 175306. 

[38] Fuchs, K. Mathematical Proceedings of the Cambridge Philosophical Society 1938, 34, 100- 


12 



108. 


[39] Sondheimer, E. Advances in Physics 1952, 1, 1-42. 

[40] Ward, A.; Broido, D. A. Phys. Rev. B 2010, 81, 085205. 

[41] Esfarjani, K.; Chen, G.; Stokes, H. T. Phys. Rev. B 2011, 84 , 085204. 

[42] Dolling, G.; Brockhouse, B. N. Phys. Rev. 1962, 128, 1120-1123. 

[43] Xiang, J.; Drzal, L. T. Carbon 2011, 49-, 773 - 778. 

[44] Veca, L. M.; Meziani, M. J.; Wang, W.; Wang, X.; Lu, F.; Zhang, P.; Lin, Y.; Fee, R.; 
Gonnell, J. W.; Sun, Y.-P. Advanced Materials 2009, 21, 2088-2092. 

[45] Tian, X.; Itkis, M. E.; Bekyarova, E. B.; Haddon, R. G. Sci. Rep. 2013, 3, 1710. 

[46] Lu, J.; Zheng, Y.; Sorkin, A.; Loh, K. P. Small 2012, 8, 3728-3732. 

[47] Dubrovkin, A. M.; Tao, J.; Ghao Yu, X.; Zheludev, N. L; Jie Wang, Q. Sci. Rep. 2015, 5, 
9837. 


13 






2 6 10 
Time (ns) 


2 6 10 
Time (ns) 



14 


























FIG. 1: Typical sample structure, experimental data and their corresponding fitting cnrves. (a) 
A three-dimensional cutaway view of a typical sample, (b) An optical image of a graphite sample 
exfoliated on top of sapphire substrate. The solid blue and red circles indicate typical sizes of pump 
and probe spots, respectively. The scale bar is 100 ^m. The (c) normalized amplitude, (d) phase, 
(e) in-phase, and (f) out-of-phase signals from the lock-in amplifier are presented as a function of 
delay time. The blue solid lines indicate the result of the fitting using the thermal model, showing 
an excellent fit. The dark and light green dash lines indicate the cases when kc is 10 % higher and 
lower than the ideal fitting result, respectively. 
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FIG. 2: Temperature dependent and thickness dependent thermal conductivities. Temperature 
dependent data (symbols connected with dotted lines) of samples thicker than 400 nm (a) and 
thinner than 400 nm (b) . Data from Ref. [12] (squares) for the same thickness range is also shown. 
Similar trends are observed. Thermal conductivity of graphite samples as a function of thickness 
for temperatures (c) greater than 100 K and (d) less than and including 100 K. Experimental data 
from Ref. [T2| at room temperature only (squares) are also shown here for comparison. Both sets 
of room temperature data present the same trend. denotes the data from Ref. |T2l 
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FIG. 3; TEM images on the cross section of a bulk graphite sample. Bright field (BF) (a) and dark 
field (DF) (b) TEM images on the same region of the same sample. Numerous wide dark strips in 
the BF image can be easily located in the DF image, which are shown as bright strips and indicate 
grains of a specific orientation distinct from the other areas, (c) A magnified BF TEM image of 
dark and bright regions. Non-uniformity within dark regions can be identified by the fluctuation 
in the contrast, (d) High resolution TEM image on the boundary of a dark and a bright region. 
Both regions demonstrate nice atomic layered structure, which indicates that they share the same 
c-axis but are of different rotational orientations. Scale bars are 1 ^m (a), 1 (b), 200 nm (c) 

and 10 nm (d), respectively. 
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FIG. 4: Thickness dependent thermal conductivities and c-MFP spectra reconstruction. Experi¬ 
mental measurements of kc (blue open squares) and the corresponding reconstructed accumulative 
thermal conductivity as a function of c-MFP (red lines) in graphite samples at various tempera¬ 
tures: (a) 294 K, (b) 100 K, (c) 60 K and (d) 40 K. As temperature decreases, the c-MFPs increase 
to a maximum of around 600 nm at 40 K. The x-axis corresponds to thickness for the experimental 
data and c-MFP for the c-MFP reconstruction. Inset: The suppression function (blue solid line) 
and the kernel function (red solid line) versus Knudsen number, which are used to perform the 
reconstruction. 
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FIG. 5: Distributions of c-MFPs in HOPG samples: (a) Differential phonon c-MFP spectra at 
different temperatures, (b) Upper (green dashed line) and lower (red dashed line) limit of phonon 
c-MFPs which carry and transport heat. The black dashed line indicates the c-MPP^id, or the 
midpoint c-MFP. The blue solid line indicates the c-MFPp, corresponding to the c-MFP where the 
maximum of the differential c-MFP spectrum occurs. The c-MFP bandwidth is much narrower 
than in isotropic crystals. 
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